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ABSTRACT 

OP I We present the Submilhmeter Array observation of the CO J=2— 1 transition to- 

f^ ' wards the northern galaxy, ARP 302N, of the early merging system, ARP 302. Our high 

CN . angular resolution observation reveals the extended spatial distribution of the molecular 

gas in ARP 302N. We find that the molecular gas has a very asymmetric distribution 
with two strong concentrations on either side of the center together with a weaker one 
QQ . offset by about 8 kpc to the north. The molecular gas distribution is also found to 

be consistent with that from the hot dust as traced by the 24 /im continuum emission 
observed by the Spitzer. The line ratio of CO J=2— 1/1— is found to vary strongly 
T^ • from about 0.7 near the galaxy center to 0.4 in the outer part of the galaxy. Excitation 

analysis suggests that the gas density is low, less than 10^ cm~^, over the entire galaxy. 
O I By fitting the SED of ARP 302N in the far infrared we obtain a dust temperature of 

Td=26-36 K and a dust mass of M(just=2-0-3.6xl0^ Mq. The spectral index of the 
C^ . radio continuum is around 0.9. The spatial distribution and spectral index of the radio 

continuum emission suggests that most of the radio continuum emission is synchrotron 
emission from the star forming regions at the nucleus and ARP302N-cm. The good 
spatial correspondance between the 3.6 cm radio continuum emission, the Spitzer 8 &: 
24 /im data and the high resolution CO J=2— 1 observation from the SMA shows that 
0^ ' there is the asymmetrical star forming activities in ARP 302N. 
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Most luminous infrared (IR) galaxies (LIRGs; 10^^ Lq < Ljr < 10^^ Lq) are found in in- 
teracting or merging systems with disturbed morphologies. In the standard scenario proposed by 
Sanders et al. (1988a, b), the formation of LIRGs, or ultra-luminous infrared galaxies (ULIRGs) is 
due to t he interaction and merging of two gas-rich galaxies. Numerical simulations of interacting 
galaxies (JBarnes &: Hernquistlll99ll : iMihos Sz Hernquistlll996l ) indicate that a stable bar is formed 
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during the first encounter. Under the gravitational torque induced by the bar, the molecular gas 
is channelled into the central regions of the galaxies f orming a huge reservoir to fuel intense star- 
bursts. Recent high angular resolution observations (JDownes &: Solomonlll998l ) seem to support 
this picture by revealing rotating disks of molecular gas that has been driven into the centers of 
the mergers. 

A lot of (U)L IRGs have bee n studied at high aii gular resol ution, for example A RP 299 



Casohet al.1ll999l l. NGC5394/9 5 (JKaufman et al 



Qlsson et al. 



Zhu et al. Ibnn^l^. N GC5194/95 (JMatsushita et al. 



20021 '! ■ ARP 220 (IWiedner et a l. 2002), NGC4038/39 
20041 ^. NGC 6090 (|Wang et al. ii2004i ). NGC5218 



2007^. A few sam 



vestigated. IXilouris et al. 



jles of the sources along the merging sequence are also been in- 
(120041 ) found evidence for an increase in 100- to 850-/im flux density 
ratio along the merging sequence and interpreted this trend as an increase of the mass fraction of 
the war m dust in comp a ,rison to the cold dust, inferring the possible enhanced star formation rate 



(SFR). iBergvall et al. I (120031). however, quest ion the ubiquitous connection between interactions 



and starburst activity. iGao &: Solomon I ( 19991 ) find the enhanced star formation efficiency (SFE) 
for the later merging systems, but the effect is cancelled by the decrease in the amount of molecular 
gas. That might result in no change in the infrared luminosity, i.e. the indicator of the starburst 
activity. Thus more sources in high resolution are needed and multiple CO transtions are essential 
to investigate the gas properties for the mergers at different stages. 

Here we present another case study of an early merger system, ARP 302. This galaxy is also 
named as VV 340 or UGC 9618, which consists of two galaxies (ARP 302N and ARP 302S) with 
a separation of about 40" between the two nuclei, which is ~27 kpc at the angular distance of 138 
Mpc (the luminosity distance is 147 Mpc), and has an infrared luminosity of 4.5x10^^ Lq. ARP 302 
is the most CO luminous system known in the local universe, and more luminous in CO than those 
ultraluminous galaxie s out to redshifts of 0.27 ( Solomon et al. 1997). The detected molecular gas 
mass is 9.1x10^*^ M© (JLo. Gao &: Gruendllll997l ) (corr ected by the adopted comsmology here). H I 



data gives a mass of atomic hydrogen 2.7x10^^^ M© (jvan Driel. Gao &: Monnier-Ragaignell200ll ) 



In all observations it was found that ARP 302N contributes most of the IR emission and molecular 

gas mass. 



Previous high resolution CO J=l— observation on ARP 302 by lLo. Gao &: Gruendll (jl997l ) 
(the angular resolution is 7"x5") reveals two extended CO disks (> 10 kpc), and ARP 302N shows 
roughly an exponential disk with a scale length of 7 kpc exte nding out to 23 kpc in diameter. Based 
on the unusually extended distribution of the molecular gas, ILo. Gao &: Gruendll ( 19971 ) concluded 
that the high IR luminosity in ARP 302 arises simply because an extremely large amount of 
molecular gas is forming high-mass stars at a rate similar to the Galaxy, instead of the starburst 
activities in the galactic nucleus. 

In this study, we report the high angular resolution CO J=2— 1 observation of ARP 302N and 
the analysis of archival data in the radio and mid/far infrared to study the physical properties of 
the gas as well as the star formation activities in this early merger. In ^ we show the details of 
our observations together with the information about the archive data used in this paper. In §3.H 
^3.21 and N3.31 we present the results of spatial distribution and the kinematics of the molecular gas 
of our observations. In §3.41 and §3.51 the results of the radio emission and the mid/far infrared 
emission from the archive data are presented. Finally we present the dust properties using a dust 



model in ^4.1^ the excitation condition of the gas in N4.21 the star formation activities in ^4.31 and 
a comparison with other mergers in ^4.4[ A summary of this study is given ii i j^ In this paper, w e 



use a A cosmology with Hq = 71 km s"^ Mpc"\ S^a = 0.73 and 0^ = 0.27 (jSpergel et al.lboosi ). 
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2. Observations 

2.1. CO J=2— 1 observation with the Submillimeter array 

We carried out the observation of the CO J=2— 1 line fr om ARP 302N on 2005 March 26, 
g the Submillimeter array (SMA^) ( Ho. Moran fc Loll2004l ) in its most compact configuration. 



At that time the SMA consists of 7 antennas of 6 meter in diameter. The weather condition during 
our observations is average, with atmospheric opacity at zenith T225 ^ 0.35. The single sideband 
system temperature is about 450 K. Because of the large extent of ARP 302N, we cover the whole 
galaxy using a two-pointing mosaic, with a separation of 20" between the two pointings located at 
10" north and south of the nucleus (a=14'*57™00^67, 6=+2A°37'02" 



1" 



The total on-source time was about 2 hours per one pointing. We observed the CO J=2— 1 line 
in the upper sideband of the double sideband receiver. The SMA digital correlator was configured to 
provide twenty-four independent spectral windows in each sideband. Each window was divided into 
128 channels, which provided 0.8 MHz frequency resolution (i.e. ^1.1 km s~^) at 223 GHz, giving 
a total frequency coverage in each sideband of 1.97 GHz (i.e. ~2650 km s^^). We tuned the SMA 
receivers to the redshifted CO J =2—1 line at the frequency of 222.967 GHz, which corresponds to 
the systemic velocity of 10180 kms~^. We observed Jupiter and Callisto for bandpass and flux 
calibration respectively, and estimated the flux density scale to be accurate to about 10% at 1.3 mm 
band. The nearby quasar 1357-1-193 was observed every 20 minutes to correct for time-dependent 
gain variations due to atmospheric fluctuations. The data were calibrated using the MIR-IDL 
package, which was specially developed for the SMA. We then exported the data into MIRIAD 
format for further processing. The mos aiced images made from SMA data were deconvolved using 



the Steer-Dewdney-Ito clean algorithm (ISteer et al. Ill984 ). The synthesized beam size is 3.4"x2.5" 



for the CO J=2— 1 observatrion with the uniformed weighting and the resulting noise level is 60 
mJy beam~^ at a velocity resolution of 20 km s~^. The 230 GHz continuum image of ARP 302N 
was formed from the line-free channels and the rms noise level is 5 mJy beam~^. No 1.3 mm 
continuum emission was detected from our observation at 3 a level of 15 mJy beam~^. 



2.2. CO J=2— 1 observation vi^ith the Submillimeter Telescope 

To estimate the amount of CO J=2-l flux recovered in our SMA observations, we also ob- 
served the CO J=2— 1 line in ARP 302N with a single dish telescope, the Submillimeter Telescope 
(SMT), in Arizona on 2005 December 14 under excellent weather condition (t225GHz ^ 0.08). The 
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observation was done in beam switching mode, with a beam throw of 2 arcmins and a chopping 
frequency of 1.25 Hz. The half-power beamwidth (HPBW) of the SMT telescope at the frequency of 
CO J=2-l line is 28", which is large enough to cover most of the CO emitting region of ARP 302N. 
The acousto-optical spectrometer (AOS) was used as the backend, with a usable bandwidth of 1 
GHz and a velocity resolution of 1.3 km s^^. A single pointing observation was made toward the 
nucleus of ARP 302N. The on-source integration time was 45 minutes. The data was reduced with 
the GILDAS software package. All the spectra were then co-added to produce the final CO J=2-l 
spectrum of ARP 302N. The final spectra were smoothed to the resolution of 14 km s^^. The 
noise level on antenna temperature scale was 3 mK. We converted the spectrum into main beam 
temperature scale using a main beam efficiency of 0.66. 



2.3. VLA archive data 

We searched the Very Large Array data archive, which is provided by the National Radio 
Astronomy Observatory (VLA®), for data on ARP 302N. Previous observations were found for X 
band (8.44 GHz, i.e. 3.6 cm), C band (4.86 GHz, i.e. 6.0 cm) and L band (1.49 GHz, i.e. 20 cm). 
We used the NRAO data reduction package to reduce the available data sets. Continuum images 
were made using a natural weighting. For X-band data, which were obtained in the extended 
configuration, we applied tappering to the data in order to recover the more extended continuum 
emission. 



2.4. Spitzer archive data 

To compare with our SMA data, we also searched for the archival data on ARP 302N produced 
by the Spitzer telescope. We found that ARP 302N was imaged with the Infrared Array Camera 
(IRAC; Fazio et al. 2004) in ah four channels (3.6, 4.5, 5.8, and 8.0 nm) on 2005 July 16 and with 
Multiband Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004) in 24, 70 and 160 //m on 
2005 January 25. The images were reduced with the Spitzer Science Center (SSC) IRAC reduction 
pipeline and combined with the SSC mosaicker. The point spread functions are 2.2", 6", 18" and 
40" for 8, 24, 70 and 160 /xm wavelengths, respectively. We only use the 8.0 /im and the MIPS 
data in this paper. The uncertainties on the final absolute calibrations are estimated at 10%, 10%, 
20%, and 20% for the 8, 24, 70, and 160 /um data, respectively. 
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3. Results 

3.1. Spatial distribution of the molecular gas 

In Figure [T]we show the channel maps of CO J=2— 1 emission for ARP 302N. There is a clear 
positional shift of the emission centroid with the velocity along the north-south direction, consistent 
with the rotation of the gas in the galactic disk. In each velocity channel we find that the emission 
is just barely resolved in the east-west direction, which is perpendicular to the galactic plane. In the 
channels with the velocity within ±80 kms~^, the emission is compact and only slightly resolved 
by the SMA, showing a symmetric distribution. The emission in this velocity range is centered at 
the nucleus and the emission of some channels is weaker than the channels outside this range. In 
higher redshifted and blueshifted velocity channels the emission is slightly more extended in the 
north-south direction, especially for the channels above ±220 kms~^. For the receding part of the 
galaxy the CO J=2-l emission is resolved into two distinct components, which are separated by 
about 8". 

Because the total extent of the CO emission is significant by comparison with the primary 
beam of the SMA, part of the CO flux could be resolved out in our SMA observations. In Figure [2] 
we compare the spectrum of CO J=2-l obtained with the SMT at the resolution of 28" and that 
obtained by convolving our SMA channel maps to the same resolution. The CO J=2— 1 flux 
detected by the SMA within one SMT beam centered at the nucleus is 356(±8) Jy km s~^. It 
suggests that the SMA recovers about 60% of the CO J=2— 1 flux from Arp 302 in comparison 
with a total flux of 576 Jy kms~^detected by the SMT. The total CO J=2— 1 flux detected by 
the SMA is 382 (±10) Jy km s~^. The most significant difference between both spectrum is the 
comparatively lower emission in the SMA spectrum from 220 to 320 kms~^, which could be due 
to extended emission that is partially resolved out by the SMA in the outer northern region. The 
similarity in the line shape between single dish and the SMA spectrum suggests that the missing 
fiux is more or less uniformly distributed across the line. 

Figure[3l^b) shows the contours of integrated CO J=2— 1 intensity map at 3.4"x2.5" resolution 
by the SMA. The molecular gas almost covers the entire optical disk of the galaxy seen almost edge 
on and distributes over a region with the size of about 15 kpc. We find that the molecular gas has a 
very asymmetric distribution with two strong concentrations on either side of center together with 
a outmost weaker concentration offset by about 8 kpc to the north. In the figure, we show three 
crosses from bottom to top: one is at the nucleus, the second is at the peak of 3.6 cm emission 
outside of the nucleus (see ^3.4p and the third one is at the outmost north for the weak CO J=2-l 
emission concentration. To compare with the BIMA data, we show the BIMA CO J=l— map in 
Figure [3^a). Our higher resolution map allows us to investigate the molecular gas distribution with 
three concentrations along the north-south direction, which the BIMA data can not tell us. We also 
convolved the SMA beam to the BIMA's and the map is shown in Figure El^d). The BIMA map 
shows more extended distribution of the gas with large radial distance comparison with the SMA 
result, which may be due to the different exciation conditions for the gas at large radial distance. 



3.2. CO J=2-l/J=l-0 line ratio 

In Figure [3||e) , we present the map of the brightness temperature ratio of CO J=2— 1 emission 
from the SMA to CO J=l— emission from the BIMA. To calculate the line ration, we first 
convolve our CO J=2— 1 map to the same angular resolution of 7"x5" of the BIMA data as shown 
in Figure md) and form the line ratio map. We find that the line ratio varies from 0.2 at the outer 
region to about 0.5 at the peaks. The ratio should be the lower limit since the missing flux of 
the SMA is about 36% but there is no missing flux for the BIMA data - the flux of 304ib45 Jy 
km s~^ from the BIMA CO J=l— observation ( Lo. Gao &: Gruendllll997l ) is similar to the value 



of 290ib35 Jy km s~^ from NRAO 12m telescope (private communication from Y. Gao; Gao & 
Solomon 1999). The line ratio is a very important observational constraint to estimate the physical 
conditions of the molecular gas in Arp 302N as we will show in the next section. 



3.3. Kinematics 

We present in Figure ^c) the intensity-weighted velocity field as contours. A velocity field 
typical of normal rotation is clearly seen along the major axis within a radius of ~4", i.e. 2.7 kpc. 
Along the minor axis, we can not detect any disturblance of the gas in the inner region. In Figured] 
we shows the position-velocity (PV) diagram of the CO J=2-l emission along the major axis 
(position angle of 1°). We note that the molecular gas exhibits a nearly constant velocity gradient 
(about 100 km s~^ kpc~^) within the central 6" or 4.0 kpc in linear scale. In the velocity range 
between 260 and 340 km s~^ there are clearly two components: one is associated with the gas 
at ARP302N-cm and the other component located to the northern part about 8" to 15" away 
from the nucleus. The gas at 3.3" (2.2 kpc) north of the nucleus shows a slightly steeper velocity 
gradient (about 170 km s~^ kpc~^) and the position of this feature is consistent with ARP302N- 
cm (the secondary 3.6 cm radio emission feature, see ^3.4p . The molecular gas located at large 
galactocentric distance toward the north shows clear deviation from solid body rotation, as can be 
seen in Figure HI 

Sakamoto et al. I ( 19991 ) find that there is a steep rise in rotation velocity, i.e. the veloc- 



ity gradient ~ 1000 km s~^ kpc~^, in the central regions with radii of typically < 0.5 kpc to- 
wards most nearby spiral galaxies. It is therefore apparent that the presence of a shallow velocity 
gradient over an extended region of 4.0 kpc in radius sets Arp 302N apart from normal cases. 
However, solid body rotatio n over kpc-s cales is found in s ome special cases, such as NGC 4194 
(JAalto fc Hiittemeisterlbood'l. NGC 7479 (iLaine et al.l ll999l'l. UGC 2866 (JHiittemeister et al.lll999l ) 



and NGC 891 (JGarcia-Burillo &: Guelinlll995l ). Usually the presence of a molecular ring or the solid 



body rotation in bars are evoked to explain the observed linear velocity gradient. 



3.4. VLA results 

Low resolution of the images at C band (6 cm) and L band (20 cm), as shown in Figure O 
indicate that the radio emission is concentrated in the central region. There is only one central 
peak and the details can not be well resolved with the current poor resolutions, which are 13-14" 



for both observations. The integrated intensities are 30.2ib0.4 mJy and 89.0±3.6 niJy for 6 cm and 
20 cm respectively, yielding a spectral index, a (where S a i^~"), of 0.9. 

High resolution of the X band (3.6 cm) images are shown in Figure [5l 3.6 cm continuum 
emission is mainly contributed from two components with similar total integrated intensity. One 
component, associated with the galactic nucleus, is compact and marginally resolved (right top 
panel of Figure E}. A second component is located at 3.3" (~ 2.2 kpc) north of the nucleus 
(ARP302N-cm hereafter) and extends from 2" to 4" away from the nucleus (left top panel of 
Figure (H) . All the results are summarized in Table [TJ 



3.5. Spitzer results 

In Figure [H we show the mid infrared continuum emision from Arp 302 at 8.0 //m and 24 /im 
taken with IRAC and MIPS instruments of Spitzer. Both 8.0 \vcq. and 24 //m images exhibit similar 
morphology to that seen in the CO J=2-l emission. Hot dust traced by the 8.0 /um and 24 /im 
share the similar spatial distribution as the molecular gas traced by the CO J=2-l emission. With 
a resolution of ~7" at 24 /xm, ARP 302 is very well resolved. From the MIPS image we estimate 
the flux densities of 0.32±0.03 Jy and 0.08±0.01 Jy for ARP 302N and ARP 302S respectively. The 
flux densities are 5.16±1.03 Jy at 70 /im and 5.08±1.02 Jy at 160 /im for ARP 302N respectively. 
Applying the same flux ratio based on 24 /tm between ARP 302N and ARP 302S to the total IRAS 
flux densities of ARP 302, we estimate 5.32ibl.60 Jy at 60 /xm and 11.6ib3.5 Jy at 100 /xm for 
ARP 302N. 



4. Discussion 

4.1. Dust properties constrained by one component dust model 

The dust continuum allows us to investigate the properties of the dust, such as the dust 
temperature and the dust mass. To model the dust emission from ARP 302N, we use the IR 
continuum fluxes measured by Spitzer MIPS at 70 /xm (5.16ibl.03 Jy) and 160 /xm (5.08ibl.02 Jy), 
and also by IRAS at 60 /xm (5.32±1.60 Jy) and 100 /xm (11.6 ifc3.5 Jy) (see detail s in gM]) together 
with the flux density of 135ib24 mJy at 850 /xm reported bv iDunne et al. I ( 200d ). 



Assuming that the dust emission is optically thin, the flux density of the dust emission is 
related to the dust temperature and the dust mass by the relation: Sjy = (1 + z)[Bi,(T(just) — 
B[y(Tbg)]Atd(z^)Mdust/DL, where S^, is the observed flux density, Tjust and Tbg are the dust and 
background temperature respectively, Bi, is the Planck function, Dl is the luminosity distance of 
the source, Mdust is the dust mass and kj is the dust absorption coefficient. For the dust absorption 
coefficient, we adopt K^iy) = 0.4 (i//250 GHz)'^ cm^ g~^, where /? is the emissivity index. We note 
that the emisivity index (5 depends on the dust properties and needs to be specified beforehand. By 
fitting the observational measurements, we can derive the dust mass and dust temperature for each 
specified value of /3. In Figure [71 we show three fitting results with (/3, Tdust, Mdust) combinations: 
/?=1.0, rdust=36 K and Mdust=2.2xl0S Mq (the best fit), /3=1.5, Tdust=29 K and Mdust=2.0xl08 



Mq (moderate f3 value), and /3=2.0, T(iust=26 K and M(just=l-lxlO^ M© (the simple modified 
blackbody fit). The dust temperature ranges from 26 to 36 K for the different models. Previous 
two models gives a dust mass of 2.0-2.2x10® Mq, which are twice of the value from the simple 
modified blackbody model. 

The simple modified blackbody model underestimates the 850 ^m flux density by ~60%, which 
indicates that either additional cold component or modifications to the emissivity index are needed 
or that the optically thin assumption is inappropriate. Actually, the gas is very extended in this 
galaxy and should result in the similar distribution for the dust emission. The dust emission is 
still optically thin at 70 /im even if the dust is distribuated over a region with a small radius of 
1 kpc, with a dust mass of 2.0x10® Mq, thus th e optically thin assu mption should be correct in 
this galaxy. In a large survey for local galaxies bv lDunne et al.l (2000), the dust emissivity indexes 



to the best fitting models are found to be 1.3ib0.2 and they argued that the galaxies may contain 
a significant amount of cold dust. For ARP 302N, an additional cold gas component with a dust 
mass of 1.6-2.5x10® Mq and a dust temperature of 20-15 K can conpensate for the lower predicted 
flux at 850 ^m. This results in a total dust mass of 2.7-3.6x10® Mq dependent on the temperature 
of the cold dust component. 



The observed flux in CO J=l-0 (|Lo. Gao k Gruendllll997l l is 304 Jy kms'^for ARP 302N 



and the corresponding CO line luminosities are 1.6x10^^^ K kms~^ pc^. This places ARP 302 to be 
the most CO luminous system known in the local universe and possibly the most gas rich system. 
Usually, we can estimate the gas mass using a standard conversion factor, 4.8 Mq (K kms~^ pc^)^^ 
in our Galaxy, and it leads to a gas mass of 7.7x10^'^Mq. We derived a dust mass of 2.0-3.6x10® 
Mq from our dust models as above. Adopting the gas-to-dust mass ratio of 150 in our Galaxy, this 
yields the value of 3.0-5.4x10^^Mq for the gas mass. Therefore, the gas mass estimates from the 
dust and the CO luminosity are comparable within their uncertainties and lead to a gas mass of 
5.4x10^^Mq (the mean value from two methods). 

The theoretical prediction of the flux density of the continuum emission at 1.3mm is less than 
42 mJy from the above dust models. The sensitivity of the continuum observation is 5 mJy beam~^, 
thus it is reasonable that there is no positive detection from the SMA observation if the emission 
is distributed over a region of over 3 beams. 



4.2. Excitation condition of the gas 

The excitation condition of molecular gas in ARP 302N can be studied using the CO J=2— 1/ J=l 
line ratio. As seen in Figure [3] the line ratio peaks at either side of the center. The ratio varies 
significantly across the galaxy, from about 0.2 in the outer region to about 0.5 at the peaks. As- 
suming that the missing fiux of the SMA observation is uniformly distributed, the line ratio will be 
0.2 higher and the ratio will vary from 0.4 in the outer region to about 0.7 at the peaks. The overall 
line ratio from the single dish measurements is about 0.5. From a survey of t he CO J=l— and 



J=2— 1 lines in the inner kpc regions of a large sample of nearby spiral galaxies, iBraine &: Combes 



(1993) measured an average line ratio of 0.89, but 7 out of 35 galaxies have CO J=2— 1/J=1— 



line ratios less than 0.6. Furthermore, the gas in the central regions tend to have high excitation 
temperature, which will result in the higher line ratio if the CO line emission is optically thick. 



The disks of NGC1808 (jAalto et al. 111994 ) and our Galaxy show values of 0.7 and 0.6, which are 
comparable to the value in the central region of ARP 302N. 

To be more quantitative, we use our molecular excitation model to constrain the physical 
conditions of the molecular gas. The transfer of radiation is calculated within the framework of the 
large velocity gradient approximation. We take into account all rotational levels of the CO molecule 
up J= 10. We use the collisional cross sections between CO and H2 molecules from lFlower &: Launav 



( 19851 ). The statistical equillibrium equations are solved using the Newton-Raphson method. In our 
molecular excitation model three parameters are required: the gas density, the kinetic temperature 
and the fractional abundance of the CO molecule per unit velocity gradient. At the expected 
low density of the molecular gas, the kinetic temperature is expected to decoupled from the dust 
temperature. In this case, the kinetic temperature might be as low as half of the dust temperature 
( Downes &: Solomonlll998l ). In the previous section, we estimate the dust temperature in the range 
of 26 to 36 K. For the excitation calculation, we adopt two representative values of 20 and 30 K. 
In Figure [8] (left panel) we present the dependence of CO J=2— 1/J=1— line ratio on the gas 
density ra(H2) and the abundance per velocity gradient [CO]/(dv/dr) for kinetic temperature of 20 
and 30 K. The figure shows that for the line ratio measured in ARP 302N the gas density is low, 
less than 10'^ cm~^. At this low density, the transition J=2— 1 is sub-thermally excited. Lower line 
ratios correspond to lower densities for a reasonable range of abundance per velocity gradient. To 
provide further support for this conclusion we also show the dependence of line ratio on the gas 
density and kinetic temperature T^in for a fixed abundance per velocity gradient of [CO]/(dv/dr) 
= Ix 10~^ pc (kms~^)~^ (right panel of Figure E]) . The result suggest s only a weak depende nce 
on kinetic temperature. Our finding is very similar to the conclusion of lSakamoto et al. I (|1994| ) in 
that the variation in line ratio mainly reflects the variation in gas density. To set more stringent 
limits on the gas density and temperature would require observations of higher transitions of CO 
and/or its ^'^CO isotope. 

The relatively low gas density in ARP 30 2N is reflected in its spatial d istribution. For ULIRGs, 
which are in a more advanced merger phase. iDownes fc SolomonI ( 19981 ) show that the molecular 
gas is more strongly concentrated in the nuclear region of the galaxy in the form of a rotating disk 
with size of 1 kpc or less. The estimated average gas density is higher, ranging from a few 10^ to 
10^ cm~^. In ARP 302N, there is a large amount of molecular gas (5.4x10^*^ Mq, see in ^4.ip . but 
this gas is distributed over a very extended galactic disk with a size of 15 kpc, resulting in lower 
average density. We note that the diffuse and very extended spatial distribution of molecular gas 
in ARP 302N is unusual among local luminous and ultralumious galaxies. However, in some cases 
such as in the interacting galaxy M51, the molecular gas traced by CO emission also covers all the 
optical disk, out to at least 12 kpc radius (jSchuster et al. 1120071 ). The origin of this distribution is 
currently unknown. 



4.3. Star forming activity 

The star formation activity can be estimated by SFE, the ratio of luminosity to molecular gas 
mass, of the galaxies. For ARP 302N, we estimate an infrared luminosity of 3.4x 10^^ Lq, separating 
the contamination of ARP 302S in the manner described in ^33.51 We can also estimate the infrared 
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luminosity from the 24 /iin data with the empirical formula L//j=10z^L,y(24//m.) ( Bell et al.ll2005l ) 



yielding a luminosity of 2.7 xlO^^ Lq, which is comparable with our estimated IRAS luminosity. 
The molecular gas estimate is 5.4x10^° Mq for ARP 302N. Its SFE of 6.3 Lq/Mq is slightly higher 
than that in our Galaxy but much lower than starburst galaxies or (U)LIRGs. The high infrared 
luminosity, the large amount of the gas mass and the low SFE make this galaxy peculiar. 

The SFR for AR P 302N is about 40 Mq yr'^ assuming SFR = 1.17x10-1° (Lfir/Lq) Mq 
yr~i ( Kennicuttlll998l ). and thus ARP 302N can retain such a luminosity over a period of at most 



2.0x10^ yrs . The system will likely evolve into advanced merger in a timescale of several tens 
million yrs ( Mihos &: Hernquistlll996l ). and may at that time become a ULIRG with a luminosity 



of 1.6 xlO^^ Lq even if it will lose half of its gas mass by then (given a SFE of 42 Lq/Mq for the 
merging systems with the separation less than 5 kpc from Gao & Solomon 1999). 



4-3.1. Constraint from the line ratio 

As discussed in previous section, our excitation analysis suggests that the low line ratio in 
ARP 302N is due to the low volume density of the molecular gas. The low gas density is entirely 
consistent with the low SFE obs erved in ARP 302N, as star formation is associated with dense gas 
rather than just total gas mass ( Gao Sz Solomonll2004l ). 



4.3.2. Radio continuum 

The thermal radio emission from HII regions is identifiable by its flat spectral index with a 
typical value of 0.1. A spectral index of 0.9 as measured for ARP 302N implied that non-thermal 
emission dominates. Synchrotron emission with a spectral index of 0.9 is likely to be dominated by 
supernova remnants with small diameters, referred to the star forming activity in ARP 302N. From 
the high resolution images at 3.6 cm, the continuum emission mainly comes from the nucleus and 
ARP302N-cm. Therefore, the nucleus and ARP302N-cm are the regions with strong star forming 
activity. 



4.3.3. Strong star forming activities at the nucleus and ARP302N-cm 

The extended infrared emission from Spitzer 8 and 24 /xm data excludes the existence of a 
strong AGN in ARP 302N and favors the star formation origin for the large infrared luminosity in 
ARP 302N. Therefore it fully supports that the star formation dominates the IR emission in the 
galgaxy. The 3.6 cm observation shows that there is strong and extended radio continuum emission 
at ARP302N-cm. The Spitzer 8 and 24 //m images peak at the same region, indicating that there 
is strong emission from the PAH and hot dust. The high resolution CO J=2— 1 image from the 
SMA reveals much molecular gas exists at ARP302N-cm. All of these results support that there 
is the asymmetrical star forming activities in the system and ARP302N-cm is the strongest active 
star forming region in the galaxy. 
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From 3.6 cm images with the different resolutions at the nuclear region, it is clear that at least 
36% of the emission is not associated with the unresolved structure in the nuclear region and 64% 
of the emission is originated from a compact region (from the VLA 3.6 cm data, the deconvolved 
size is 0.23" xO. 09", i.e. 143 pc x 56 pc). As a whole, the star formation is not active in ARP 302N, 
but there is strong star forming regions at the nucleus and ARP302N-cm. 



From the previous optical spectroscopy observation (IVeilleux et al. I Il995l ) , the system was 
classified as LINER type galaxy. lAlonso-Herrero et al.l (2000) claim that LINER-like optical spectra 
would be the consequence of the shock-heating by supernovae of t he surrounding gas after a burst 
of star formation, as observed in NGC 5218 ( Olsson et al.l 120071 ). This scenario to explain the 
LINER type spectrum is consistent with our radio observation which suggests that there are many 
supernova explosions. Of course, we can not exclude the possible existence of a weak AGN in the 
center. 



4.4. Comparison with other mergers 

In Arp 302N the distribution and kinematics of molecular gas show no discernable disturbance. 
That suggests the Arp 302 system is in the early merging stage. Recently, many authors studied 
the molecular gas in high angular resolution to investigate the effects of the merging procedure 
and the star formation therein. The molecular gas is often concentrated on the overlap region or 
at least have anot her peak except for the centers o f the systems with the well seper ated nuclei 
such as ARP 299 dCasoH et al. Ill999l ). NGC 6090 dWang et al. II2004I ) and VV 114 dlono et al. 
20041 ). In these cases, the starburst acitivites is prominent and the SEE is much high er than that 



i n our Galaxies. NGC 6670 is another widely seperated merger and is well studied by lWang et al. 



(J200ll ). They also found that the stars and the molecular gas do not show any clear sign of the 
merging process. However, the atomic gas is found to concentrate in the region in between the 
two galaxies, a very clear signature of the merging process. The atomic gas will be helpful to 
determine the merging process in ARP 302. Indeed, simulations show that the high SERs and the 
SEEs arise only in the late stages of the merger and the star formation does not deplete much gas 
before it reaches the central regions in the mergering systems iMihos &: Hernauist ( 1994 . 19961 ). It 
is consistent with the observed enhanced SEE in the late mergers by lGao &: Solomon I ( 19991 ). which 
suggests that the fact before its final merging stage is the reason why the SEE in this LIRG is just 
comparable to the value in our Galaxy. The huge infrared luminosity in ARP 302N is just due to 
its large amount of the gas mass instead of the starburst activity. 



5. Conclusion 



In this paper we present the SMA observations of the CO J=2— 1 transition towards the 
northern galaxy, ARP 302N, of the early merging system, ARP 302, to study the physical properties 
of the gas. To well understand the star formation activities, the VLA data at 3.6 cm, 6.0 cm, 20 cm 
and the Spitzer (both IRAC and MIPS data) archive data are also shown here. Our main results 
are summarized as the following: 
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1. Our high angular resolution CO J=2— 1 observation reveals the extended spatial distribution 
of the molecular gas in ARP 302N. We find that the molecular gas has a very asymmetric 
distribution with two strong concentrations on either side of the center together with a weaker 
one offset by about 8 kpc to the north. This spatial distribution is very similar to the 
morphology seen previously in the CO J=2-l line. 

2. The distribution of the molecular gas is also found to be consistent with that for the PAH 
and the hot dust as traced by the 8 nm and 24 ^m continuum emission observed by Spitzer. 

3. The line ratio CO J=2-l/l-0 shows that the line ratio varies strongly from about 0.7 near 
the galaxy center to 0.4 in the outer part of the galaxy. Excitation analysis suggests that the 
gas density is low, less than 10^ cm^^, over the entire galaxy. 

4. By fitting the far infrared part of the SED of ARP 302N, we obtain a dust temperature T^ 
in the range 26-36 K and a dust mass of M(iust=2-0-3.6xl0® Mq. 

5. We determine the spectral index of the radio continuum to be approximately 0.9. The spatial 
distribution and spectral index of the radio continuum emission suggests that most of the 
radio continuum emission is synchrotron emission from the star forming regions at the nucleus 
and ARP302N-cm. The good spatial correspondance between the 3.6 cm radio continuum 
emission, the Spitzer 8 and 24 /im data and the high resolution CO J=2— 1 observation from 
the SMA shows that there is the asymmetrical star forming activities in ARP 302N and 
ARP302N-cm is the strongest active star forming region in the galaxy. 



We would like to acknowledge the anonymous referee for the valuable discussions and com- 
ments. We also thank the SMA staff for maintaining the operation of the array. The SMT is 
operated by the Arizona Radio Observatory (ARO), Steward Observatory, University of Arizona. 
Y.-P. Ao thanks the ASIAA for supporting his stay as a Visiting Scholar, when much of this work 
was done. Y.-P. Ao also acknowledge the support of NSFC grant 10733030. 



REFERENCES 

Aalto, S., Booth, R. S., Black, J. H., Koribalski, B. & Wielebinski, R., 1994, A&A, 286, 365 

Aalto, S. & Hiittemeister, S., 2000, A&A, 362, 42 

Alonso-Herrero, A., Rieke, M. J., Rieke, G. H. & Shields, J. C, 2000, ApJ, 530, 688 

Barnes, J. E. & Hernquist, L. E., 1991, ApJ, 370, L65 

Bell, E. F., et al. 2005, ApJ, 625, 23 

BergvaU, N., Laurikainen, E. & Aalto, S., 2003, A&A, 405, 31 

Braine, J., k Combes, F., 1992, A&A 264, 433 



13 



Casoli, F., Willaime, M.-C, Viallefond, F. & Gerin, M., 1999, A&A, 346, 663 

Downes, D. & Solomon, P. M., 1998, ApJ, 507, 615 

Dunne, L., Eales, S., Edmunds, M., Ivison, R., Alexander, P. & Clements, D. L., 2000, MNRAS, 
315, 115 

Fazio, G. G., et al., 2004, ApJS, 154, 10 

Flower, D.R., Launay, J.M., 1985, MNRAS, 214, 271 

Gao, Y. & Solomon, M., 1999, ApJ, 512, L99 

Gao, Y. & Solomon, P. M., 2004, ApJ, 606, 271 

Garcia-Burillo, S. & Guelin, M., 1995, A&A, 299, 657 

Ho, P. T. P., Moran, J. M. & Lo, K. Y., 2004, ApJ, 616, LI 

Hiittemeister, S., Aalto, S. & Wall, W. F., 1999, A&A, 346, 45 

lono, D., Ho, P. T. P., Yun, M. S., Matsushita, S., Peck, A. B. & Sakamoto, K., 2004, ApJ, 616, 
L63 

Kaufman, M., Sheth, K., Struck, C. et al. 2002, AJ, 123, 702 

Kennicutt, R. C., 1998, ARA&A, 36, 189 

Laine, S., Kenney, J. D. P., Yun, M. S. & Gottesman, S. T., 1999, ApJ, 511, L709 

Lo, K. Y., Gao, Y. & Gruendl, R. A., 1997, ApJ, 475, L103 

Matsushita, S., Sakamoto, K., Kuo, Cheng- Yu et al. 2004, ApJ, 616, L55 

Mihos, J. C. k Hernquist, L. E., 1994, ApJ, 425, L13 

Mihos, J. C. & Hernquist, L. E., 1996, ApJ464, 641 

Olsson, E., Aalto, S., Thomasson, M., Beswick, R. & Httemeister, S., 2007, A&A, 473, 389 

Sakamoto, S., Hayashi, M., Hasegawa, T., Handa, T., Sz Oka, T., 1994, ApJ 425, 641 

Sakamoto, K., Okumura, S. K., Ishizuki, S. & Scoville, N. Z., 1999, ApJ, 525, 691 

Sanders, D. B., Soifer, B. T., Elias, J. H., Madore, B. F., Matthews, K., Neugebauer, G. &: Scoville, 
N. Z., 1988, ApJ, 325, 74 

Sanders, D. B., Soifer, B. T., Elias, J. H., Neugebauer, G. & Matthews K., 1988b, ApJ, 328, L35 

Schuster, K. F., Kramer, C, Hitschfeld, M., Garcia-Burillo, S. & Mookerjea, B., 2007, A&A, 461, 
143 

Solomon, P. M., Downes, D., Radford, S. J. E. & Barrett, J. W., 1997, ApJ, 478, 144 



14 



Spergel, D. N., Verde, L., Peiris, H., etal. 2003, ApJS, 148, 175 

Steer, D. G., Dewdney, P. E. & Ito, M. R., 1984, A&A, 137, 159 

van Driel, W., Gao, Y. & Monnier-Ragaigne, D., 2001, A&A, 368, 64 

Veilleux, S., Kim, D.-C, Sanders, D. B., Mazzarella, J. M. & Soifer, B. T., 1995, ApJS, 98, 171 

Wang, J., Zhang, Q., Wang, Z., Ho, P. T. P., Fazio, G. G. & Wu, Y., 2004, ApJ, 616, L67 

Wang, W., Lo, K. Y., Gao, Y. & Gruendl, R. A., 2001, AJ, 122, 140 

Wiedner, M.C., Wilson, CD., Harrison, A., Hills, R. E., Lay, O. P. & Carlstrom, J. E., 2002, ApJ, 
581, 229 

Xilouris, E. M., Georgakakis, A. E., Misiriotis, A. & Charmandaris, V., 2004, MNRAS, 355, 57 

Zhu, M., Seaquist, E. R. & Kuno, N., 2003, ApJ, 588, 243 



This preprint was prepared with the AAS lATpjX macros v5.2. 



15 



Table 1. Basic parameters for the SMA and VLA observations towards ARP 302N 



Data 


Freq. 


Line/cont. 


beam size, P.A. 


Image rms 


Intensity*^ 




GHz 




arcsec, degree 


mJy beam~^ 


mJy^ 


SMA 


222.967 


CO J=2-l 


3.4x2.5, 154° 


60^ 


382±10 


VLAL 


1.49 


continuum 


14x13, 140° 


1.6 


89.0±3.6 


VLA C 


4.86 


continuum 


15x13, 178° 


0.21 


30.2±0.4 


VLA X (taper 150 kA) 


8.44 


continuum 


0.95x0.90, 59° 


0.14 


1.70±0.31;2.53±0.48 


VLA X (no taper) 


8.44 


continuum 


0.24x0.23, 87° 


0.05 


1.08±0.12;0.93±0.22 



'^For VLA 3.6 cm data, the first column for the nucleus and the second column for ARP 304N-cm. 
For the SMA result, the integrated intensity is in unit of Jy km s^^. 
^Velocity resolution of 20 km s~^. 
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Fig. 1.— Channel maps of CO J=2-l line of ARP 302N. The contours are -3, 3, 4, 5, 6, 7 and 
8 x 60 mJy beam~^ (1 sigma level), with a synthesized beam of 3.4"x2.5". Three cross symbols 
from bottom to top show the positions of nuclear region, ARP302N-cm (see §3.4p and outmost 
northern CO J=2— 1 emission peak respectively. The map is centered on the nucleus of ARP 302N 
(a(J2000)=14^^57"'00^67, (5(J2000)=+24°37'02".8). 
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Fig. 2. — (a). At the position of the nucleus of ARP 302N, spectrum of CO J=2— 1 hue obtained 
with the Sub-milhrneter Telescope at a resolution of 28". The measurement is on the main beam 
brightness temperature, (b). Spectrum of CO J=2— 1 line at the same position obtained from the 
SMA convolved to the same angular resolution as the single dish observation. 
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(o) BIMA CO(l-O) 



(b) SMA C0(2-1) (c) CO velocity field (d) Convolved C0(2-1) (e) C0(2-l)/(1-0) 
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Fig. 3.— (a). Integrated intensity map of CO J=l-0 line of ARP 302N. CO data is from BIMA 
observation, obtained by Lo et al 1997. The contours shown are 6, 9, 12, 15, 20, 30, 50, 70, 90, 110, 
130, 150 and 170 x 0.6 Jy km s"^ (b). Integrated intensity map of CO J=2-l hne of ARP 302N. 
The contours shown are 3, 6, 9, 12, 15, 20, 30, 50, 70, 90, 110 and 130 x 1 Jy km s^^ (c). CO 
J=2— 1 velocity field obtained by the SMA. Contour levels are from -250 to 320 kms~^ by the 
step of 30 kms~^ . Dashed line is blueshifted velocity and solid line redshifted. (d). CO J=2— 1 
integrated intensity map convolved to the BIMA resolution. The contours shown are 3, 6, 9, 12, 
15, 20, 30, 50, 70 and 90 x 2 Jy km s'^ (e). Ratio map of rcoj=2— i/rcoj=i— o of ARP 302N. 
The contour is from 0.2 to 0.45 by the step of 0.05. The cross symbols are the same as in Figured] 
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Fig. 4. — Position-velocity diagram of ARP 302N along major axis. The contour levels are from 
0.12 to 0.54 Jy beam~^ by a step of 0.06 Jy beam~^. 
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Fig. 5. — Radio images are shown towards ARP 302. The images at 3.6 cm are shown at two 
top panels with the resolutions of 0.76x0.54 (left top panel) and 0.23x0.09 (right top panel). The 
images at 20 cm and 6 cm are shown at bottom panels with resolutions of 13.7" x 12.7" (left bottom 
panel) and 14.8" x 13.3" (right bottom panel). The noise values and contour levels are shown at the 
right bottom of each panel. The cross symbols are the same as in Figure [TJ 
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Fig. 6. — The Spitzer 8 /im and 24 //m images are shown in contours and overlaid on the SMA CO 
J=2— 1 integrated intensity in grey scale. The contour levels are from 5% to 95% by a step of 10% 
of the peak values. The cross symbols are the same as in Figure [H 
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Fig. 7. — Observed SED of ARP 302N fitted by different single component dust models. The single 
component dust model are shown for (/?, Tjust; Mjust) combinations as follows: the best fit (solid 
line: 1.0, 36 K, 2.2x10® Mq); moderate (3 value (dashed line: 1.5, 29 K, 2.0x10*^ Mq); the simple 
modified blackbody fit (dotted line: 2.0, 26 K, 1.1x10® Mq). The vetical line shows the predicted 
value range at 1.3 mm. Note that 24 ^um data are not used for the fitting. 
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Fig. 8. — The dependence of CO J=2— 1/J=1— line ratio on the gas density n(H2) and the 
abundance per velocity gradient [CO]/(dv/dr) for the kinetic temperature of 20 K and 30 K (two 
left panels) , and the dependence of line ratio on the gas density and kinetic temperature Tkin for a 
fixed abundance per velocity gradient of [CO]/(dv/dr) = Ix 10^^ pc (kms^^)^^ (right panel). Te 
solid lines show various values of the line ratio between 0.4 and 0.8 by step of 0.1. 



